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Abstract. — We show how the interface curvature autocorrelation function 
(ICAF) and associated structure factor (ICSF), of relevance in non-equilibrium 
pattern- formation problems where sharp interfaces are present, provide new and 
interesting information on domain structure, as yet not visible via the order- 
parameter structure factor (OPSF). This is done by discussing numerical simu- 
lations of model A (non-conserved relaxational phase-ordering kinetics) in two- 
dimensional systems. The ICAF is Gaussian over short distances and exhibits 
dynamical scaling and t 1 / 2 power-law growth. We use it to show what the typical 
length-scale in the model A dynamics corresponds to physically and how it can 
be obtained uniquely, rather than simply within a multiplicative constant. Ex- 
perimental methods to measure the ICAF and/or ICSF are still needed at this 
point. 



1 Introduction 

Dynamics of pattern-formation in non-equilibrium systems is a very challenging prob- 
lem HJ . It is ubiquitous in nature and its understanding is of interest and importance 
in physical and biological sciences. It is customary to use light or particle scattering 
experiments on such systems to measure the spatial or temporal correlations, yielding 
a scattering function or structure factor. In systems which can be described by an 
order-parameter, the order-parameter structure factor (OPSF) has traditionally been 
the quantity of interest, as it is easily measurable experimentally and is well-defined 
mathematically, making comparison between theory and experiment possible. In prin- 
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ciple, the OPSF contains much of the structural information of the system's state at 
a given time. However, simulations of the well-known Q model A system, described 
below, show that some fundamental configurational information is not readily (if at 
all) visible with the OPSF. Experimental methods for measuring the ICAF or ICSF 
do not seem to exist at present. Results discussed here show a clear need for such 
experimental measurement methods. 

Hohenberg and Halperin in the late 70's proposed a classification for several types 
of pattern- forming dynamics for which a field-theoretic description existed [gj. One 
of the classes, labeled model A, is that of dissipative dynamics for a single uncoupled 
non-conserved order-parameter <j). The order-parameter could be for instance the lo- 
cal magnetization in an Ising ferromagnet, i.e. an idealized ferromagnet in which the 
magnetization of any molecule can take only two values, either +1 (up) or —1 (down). 
Model A is described by a Time-Dependent Ginzburg-Landau (TDGL) equation re- 
lating the temporal and spatial variations of the order-parameter of the system: 

M-^ = 0-0 3 + fV 2 (1) 

where M and £ are positive phenomenological constants determining the time-scale 
and interface length-scale of the dynamics, respectively. Some phenomenological pa- 
rameters have been scaled out. In this equation, there is no constraint on the average 
order-parameter per unit area as a function of time. This differs from the well known 
model B which describes spinodal decomposition in binary mixtures Q , and for which 
the order parameter is conserved. 

Under appropriate conditions, such as a critical quench, model A dynamics is 
characterized by the formation of convoluted, interpenetrating domains of two phases. 
The domains are separated from one another by sharp interfaces, i.e. sharp on the 
length-scale of the domains but smooth and of finite-width on the length-scale of the 
molecules of the system ||. Once the interfaces have formed, the system enters the so- 
called scaling (i.e. late stage) regime, where the dynamics is strongly non-linear. Note 
that the width of interfaces remains approximately constant throughout the late-stage 
regime, and is roughly 5£. In this regime, the system seeks to decrease the amount of 
interface via interface motion. 

There are several interesting aspects to the dynamics during the late-stage regime, 
but for the purposes of this article we recall only two of them. First, experiments 
and numerical simulations observe self-similar dynamics, whereby the system state 
at a given time is statistically the same as that at a later time, if space is properly 
rescaled. The minimal condition for this is that the dynamics must have the same 
time-dependence on all length-scales, so that all dynamical lengths can be expressed 
in terms of one arbitrarily chosen reference length-scale L{t) which encompasses the 
unique time-dependence. The second important characteristic of model A dynamics 
is that L(t) ~ t 1 / 2 . 
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Mathematically, the OPSF is written 

St$,t) = j($(k,t)$(-k,t)) (2) 

where A is the system area, (f>(k, t) is the Fourier transform of the order-parameter 
cp(x, t) at time t, and the angle brackets denote, as usual, an ensemble average over 
all possible initial system configurations. SAk, t) gives the statistical intensity of each 
mode of the order parameter field eft, i.e. the importance of domains whose coarse size 
is 2ir/\k\. S<f,(k,t) itself is the Fourier transform of the two-point equal-time order- 
parameter correlation function 

G<j,(x,t) = ((j)(xo + x,t)(f)(xo,t)) , (3) 

which is a measure of the spatial correlations in the order-parameter field at a given 
time. The main characteristics of S^(k,t) for model A are a maximum at k = 0, 
and for large k, ~ where d = dimensionality of space. This power-law 

tail is known as Porod's law and is a direct consequence, as shown by Porod 0, of 
the sharpness of interfaces. The k = dominant mode indicates that the probability 
of finding larger domains increases with the (coarse) size of the domain, and hence 
there is no typical coarse size for the domains. This is in stark contrast to model B 
dynamics, where S^k,t) shows a peak at non-zero k. 

The most important consequence of self-similar dynamics is that Sj,(k,t) taken at 
different times during the scaling regime can be rescaled in amplitude and k, using 
L(t) described earlier, to fall on one, universal and time-independent curve S^(y), with 
y = \k\L{t) (hence the denomination of "scaling" regime). L(t) can be computed, for 
instance, from any moment of S^,(k,t). Regardless of how it is computed, it will 
always have the same time dependence, within a multiplicative factor. The unique 
time dependence has lead to the conclusion that a typical, or dominant, length scale 
exists in model A dynamics. But the arbitrariness of the multiplicative factor has 
eluded attempts at computing or even identifying its physical origin from the order 
parameter correlation functions. Here, we solve this problem by focusing on the 
dynamics at the domain interfaces. 

Indeed, the late time dynamics of model A is known to be dominated by the 
motion of interfaces which decouple from the interior of the domains. The interface 
dynamics is curvature driven according to the Allen-Cahn result ||: the interface 
velocity at a given interfacial point is proportional to its curvature at that point; the 
proportionality coefficient is M £ 2 where the phenomenological constants M and £ 
are respectively the mobility and interfacial width (assumed time independent at late 
times). While much has been written about this interface dynamics, it pertains only to 
the behaviour normal to the interface. In contrast, correlations along the interfaces as 
they evolve have to our knowledge not been explored. We report here for the first time 
some new results for the space and time dependent curvature-curvature correlations 
along the interface. 
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Figure 1: Parametric representation of an interface in the plane, as a function of 
arclength s along the interface. 

2 Results 

An interface I in the plane can be represented parametrically by a set of vectors 
R(s) = (x(s),y(s)) defined relative to an arbitrary origin, as in fig. |l[ The parameter 
s is the arclength along the interface. In this notation, the curvature at point s is 

as^ 

The generic definition of the ICAF is 

G K (s,t) = (K(0,t) ■ K(s,t)) (4) 

where K(s,t) is the curvature of an interface as a function of arclength position s 
along the interface and time t. Numerically there are at least two ways of computing 
G K (s,t), assuming the interfaces are discretized homogeneously with a mesh of size 
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As. They both use eq. (^|) for interface i, which takes the form 



G Ki (s,t) = — - Ki(s 3 ,t) ■ Ki(s 3 + s,t). (5) 
n ^ > j=i 

where s is assumed an integer multiple of As, the subscript i refers to interface i, nj(i) 
is the number of points on the interface, and Ki(sj,t) is the curvature of the interface 
at some point Sj = jAs, at time t. 

The first method is a simple weighted average of all G K .(s,t), giving more impor- 
tance to longer interfaces. We denote it by G^(s,t), 

Gl{s,t) = Y! ^(^M) ( 6 ) 

£i=i i=i 

where the prime superscript in the sum indicates that only interfaces longer than a 
length of 2s are used and N[{t) is the number of interfaces used. 

The second definition, which uses the same notation as G^(s,t), is 

i.e. G^(s, t) is a weighted average of the normalized ICAF for each individual interface. 
Hence G^(s,t) < 1 for all s, with exact equality at s = 0, and can be interpreted as 
the average relative value of curvature a distance s on either side of a point where the 
curvature is K. Equation (|6|) turned out to be more suitable for analytical calculations, 
while eq. (Q) gives smaller statistical error numerically, but both measure the same 
correlation of local interfacial curvature. 

Numerical simulations of model A were done for flat Euclidean systems of sizes 
100 x 100 and 200 x 200, with periodic boundary conditions. Forty random order- 
parameter configurations were generated and evolved by integrating eq. (|l|) using a 
standard Forward- in-time/Centered-in-Space Euler integration scheme Q, from t = 
+ to t = 1000. The system mesh size Ax used was 1, and the time step At = 0.03. 
The results were checked to be independent of system-size. Average computation time 
required was 8 hours for the 40 runs on an HP735. A typical interface configuration 
is shown in fig. ||. 

Computing K is very difficult if the bulk description (eq. (|l|)) is used, as it requires 
extracting interfaces from the order-parameter configurations by systematically scan- 
ning these and finding all interfaces, splining them for smoothness and finally comput- 
ing K at regular intervals along an interface. An interface description, which evolves 
the interfaces directly via an interface equation 0, allows for direct computation of 
K and therefore G K (s,t), while the runtime can be decreased by a factor of 10 for 
flat systems and 50 to 100 for the curved (i.e. non-Euclidean) model A ]?]]. Such a 



5 




Figure 2: Domain interfaces in a typical two-dimensional model A system. 



discretized interface description produced results closer to the analytic predictions 
than did the discretized bulk description (c/. Discussion section), but both yielded 
otherwise identical results. 

The ICAF was computed at several different times during the scaling regime, 
t = 17 to t = 300. The rescaled ICAF, denoted 

^(xjH^^r 1 *,*), x = L (t)-\ (8) 

is shown in fig. [3L where the scaling length-scale L a (t) was arbitrarily defined as the 
first zero of G^.(s, t). The vertical error bars, not shown for clarity, are at the origin 
and increase roughly linearly to 0.01 in the vicinity of the minimum, then further 
increase to 0.02 at s = 200. The error was computed by making an analogy between 
the curvature K and magnetization rh of one-dimensional Ising magnets of different 
lengths ||. The error for distance s is then the weighted average of the deviation of 
each magnet's value of (m(O)m(s)) from the value of (m(O)m(s)) for the ensemble of 
magnets. G^(s,t) corresponds to (m(O)m(s)). This was deemed the most reasonable 
method of error calculation, given the values of K along an interface, and therefore 
the statistical error in products of K, are correlated. 
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Figure 3: E^(x), from bulk configurations at various times. L a (t) was defined as the 
first zero of G^(s,t). 

3 Discussion 

The salient features of the curves in fig. |3| are the nice superposition of the 
(within error bars), the power law for L Q (t), the negative autocorrelation for s > L (t), 
and finally the Gaussian form for distances much smaller than L (t). The perfect 
dynamical scaling indicates that the ICAF correctly captures this very important 
characteristic of model A dynamics, and that the number of runs and system size 
used give an accurate and reliable measure of G^(s,t). The power law in L D (t) is 
found through linear regression to be 0.45 ± 0.02, very close to the theoretical value 
H of 1/2 (the interface description gives a power-law even closer to the analytic 
prediciton of 1/2, namely 0.48 ± 0.01, as it is less sensitive to discretization effects). 
G^-(s, t) therefore also captures to a high accuracy the well-known power-law behavior 
characteristic of model A dynamics. 

Over short distances (up to roughly L Q (t)/2), G K (s,t) can be checked to be Gaus- 
sian. A Gaussian form for a correlation function was shown by Porod to be a charac- 
teristic of fluctuating systems without sharp variations A plot of K(s) (not shown) 
indeed looks very much like a snapshot of a one-dimensional fluctuating "membrane" , 
smooth and without any sharp variations (cf. fig. §). A first attempt at obtaining 
an approximate analytical expression for G^(s,t) is exposed here, based on a method 
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developed by T.M. Rogers in his Ph.D. thesis |10[ , in the context of model B, but 
never tested numerically. 

Consider the model A curvature equations in a flat system: 



92K + K* (9) 



at J a ds 2 

= -VaK 2 (10) 



where s is the arclength along an interface, a is a parameterization of the interface 
in which every point has constant coordinate a (exploiting the fact that the interface 
moves perpendicularly to itself [J7j]), and g is the metric on the interface, relating the 
elements of length in both gauges (s and a): 

ds = yfg da. (11) 

Now let us make the following mean-field approximation: 



~ -y/gh{t) (12) 



h{t) = ^Jx 2 ds (13) 

where L is the total length of the interface. This approximation becomes exact for 
circular domains. Furthermore, neglect the cubic term of eq. (||). This term is domi- 
nant for circular domains. For convoluted domains, numerical testing indicates K 3 is 
comparable to the diffusion term for as many as half the interface points. Therefore, 
the two approximations work in opposite directions, one becoming exact for circular 
domains, the other getting better for convoluted domains. On short length-scales, 
however, convoluted domains are locally circular, but the K 3 term should be neg- 
ligible since convoluted domains see their curvature decrease rather than increase. 
Equation (P) then becomes 



dK\ 1 d 2 K 



(14) 



8t J a g(t) 8a 2 ■ 

Going to Fourier space and making use of a change of variable for time, g(t)dt' = dt, 
the integration can now be performed, and the curvature structure factor Xq obtained: 

X q (t>) i < K q (f)K_ q {t') >= <K <§ )2> e-*V (15) 

where N is the number of points on the interface and q is the wavenumber in the 
reciprocal space of a. Assuming all K q have equal amplitude at t = 0, a backwards 
Fourier transform yields 

G 1 K (a,t') = Me-« 2 W. (16) 
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Now t' must be found as a function of t. This can be done by noting that h{t) is equal 
to G} ( (0,t'(t)), i.e. h(t) = y^. Also, 



t'(t) - 

The equation for the metric is therefore 



dt 



-1/2 



This has g(t) = t~ 1 /ir as solution, so that t'(t) = irt 2 /2. From eq. (11), a = s/^/g 



*l 2g rALw) ] (17) 



syf-nt. Substituting in eq. fllq), 



G 1 K {s,t) = ± t e- s2 "\ (18) 



which is Gaussian, has time-dependencies consistent with power-law growth in model 
A dynamics (the amplitude of G^(s,t) has units length squared), and dynamically 
scales, but does not capture the negative autocorrelations at longer distances. Quan- 
titatively, numerical simulations find the amplitude of the Gaussian to go as [2(1.1 db 
0.1)t] _1 , in very good agreement with eq. (18), whereas the width is smaller than the 



prediction (|lg) by about a factor of two. 

Of course, the two approximations that allowed for the calculation of G K (s,t) are 
reasonable only on short length-scales, so it is not surprising that eq. ([18]) does not 
capture the dominant wavelength of undulations of the interfaces. Also, all K q (0) were 
assumed equal. This is obviously wrong, since even at the earliest times the curvature 
structure factor (the Fourier transform of the G K (s, t)) shows a well-defined peak at a 
non-zero q mode. If the correct X«(0) is used, then the analytical G K (s,t) will have a 
dip at least at early times. Therefore the strongest approximation may yet lie in the 
X<j(0) rather than the mean-field and linearization approximations, though this seems 
unlikely. 

The most interesting feature of the ICAF is undoubtedly the relatively large nega- 
tive autocorrelation apparent at distances s > L (t). Hence the rescaled ICSF Sk(v) 
(Fourier transform of eq. (]|)) shows a well-defined peak at a non-zero value of wavevec- 
tor y = 2irL (t)/s, as seen in fig. The functional form can be fitted very nicely to 
ay b e~ cy for < y < 5, with a = 110 ± 3, b = 3/2 ± 0.01 and c = 4/3 ± 0.04, though 
no theoretical justification for this is known to us at present. The maximum of this 
function occurs at y = b/c = 9/8 ± 0.04, which is extremely close to 1. The null value 
of 8^(0) stems from the null area under the G^(s, t) curves, itself a direct consequence 
of 

'k { ds = (19) 



for any interface i (when the system has periodic boundary conditions), a constraint 
akin to a conservation law. The peak in §k(u) indicates that there is a dominant 
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Figure 4: The rescaled curvature structure factors §i<-(y), as a function of y = 
2-irL(t)/s, for the same times as in fig. [3| (joined into one curve for clarity). Fit 
gives a = 110, b = 3/2 and c = 4/3 with 3% error. 

wavelength for undulations of the interfaces. This dominant spatial mode is a feature 
which remains as yet either hidden in or inaccessible via the order-parameter structure 
factor. To our knowledge, it is the first time that it is clearly identified. The closeness 
of y to 1 also indicates that the first zero of G^(s,t) corrresponds to the dominant 
spatial oscillation mode for curvature undulations, i.e. the typical arclength distance 
from a point on the interface at which the curvature changes sign. 

The dominant length-scale L Q (t) for model A dynamics thus seems to have a 
different nature than, for instance, that of model B dynamics, whose OPSF itself 
shows a peak at non-zero wavevector. In the literature on model A and B dynamics 
one loosely speaks of this dominant length-scale as an average or typical domain 
size. However the concept of domain size is well defined only when the domains are 
morphologically disconnected or, if not, if they have a well-defined width. As discussed 
earlier in this article, model B satisfies this, but for model A only sufficiently off-critical 
quenches create domains whose bubble morphology lends itself to the definition of a 
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Figure 5: A domain in model A. The width of the OPSF should give an indication of 
how often domains of size A occur. The line B is what the dominant length-scale is, 
as given by the ICAF: roughly, the dominant radius of curvature. 



"typical domain size". The ICAF for model A shows that a dominant length-scale 
is present in the spatial undulations of the interfaces rather than in the size of the 
domains. The difference is schematized in fig. ||. We have referred to this dominant 
length as L (t), to differentiate it from the degenerate L(t). An interesting question 
this raises is whether model B interfaces would exhibit not only the dominant length 
of domain size, but a dominant length of interface undulation as well. 

The presence of a dominant undulation wavelength implies that model A interfaces 
may be called "random" only approximately p| , since truly random interfaces have 
an uncorrelated ICAF rather than the one found here. Several analytical methods 
developed to derive the scaling function for the OPSF make use of Gaussian assump- 
tions about the order-parameter field as well as the randomness of the interfaces [11]. 
Though for model A these gaussian assumptions appear to work, the non-gaussian cur- 
vature correlations may provide a clue into the break-down of gaussian assumptions 
for the more complex model B. 

An important difference between the OPSF and the ICSF is that the latter distin- 
guishes between the two domain morphologies of model A: bicontinuous for a critical 
quench, bubble for strongly off-critical quenches. Indeed, the OPSF is qualitatively 
the same for both morphologies, whereas for the bubble morphology the negative cur- 
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vature autocorrelations in the ICAF are non-existent: the peak in the ICSF shifts to 
y = for sufficiently off-critical quenches. The difference may be due to the absence 
of phase information in the OPSF. 

4 Conclusion 

We discussed several features of model A interfacial dynamics via the ICAF (interface 
curvature autocorrelation function) G K (s, t). Two notable characteristics of the ICAF 
were the gaussian form near s = 0, and the dip below zero beyond a certain distance. 
Though time-dependent lengths can be defined from the OPSF, and the existence 
of a typical and uniquely defined dynamical length L Q {t) can be infered from the 
universality of the OPSF, the peak at zero wavenumber prevents us from uniquely 
deducing L (t) and identifying its physical meaning. However the ICAF, through the 
existence of a zero, or equivalently the ICSF, through the location of its maximum, 
clearly answers both questions. The scaling length-scale for model A systems is hence 
in the interface undulations of the domains. This suggests that the scaling length- 
scale in model A dynamics is of a different nature than the domain-size length-scale 
of model B dynamics, where the order-parameter is conserved. It is also clear that 
the OPSF hides some interesting characteristics of model A interface dynamics and 
its domain morphology. Experimental methods for measuring the interface curvature 
autocorrelations could become useful in better characterizing the dynamics of some 
pattern forming systems where sharp interfaces are present. 
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